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Abstract. The heat equation is considered in the complex system consisting of many small 
bodies (particles) embedded in a given material. On the surfaces of the small bodies a Newton- 
type boundary condition is imposed. An equation for the limiting field is derived when the 
characteristic size a of the small bodies tends to zero, their total number Af(a) tends to infinity 
at a suitable rate, and the distance d = d(a) between neighboring small bodies tends to zero 
a « d. No periodicity is assumed about the distribution of the small bodies. 
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1. Introduction 

Let many small bodies (particles) T> m , 1 < m < M, be distributed in a bounded domain T> C R 3 , 
diaml? m = 2a. The small bodies arc distributed according to the law 

N(A) = ^—[ N(x)dx[l + o(l)],a^0. (1.1) 



Here A C T> is an arbitrary open subdomain of T>. k 6 (0,1) is a constant, N(x) > is a continuous 
function, and Af{A) is the number of the small bodies V m in A. Let us assume that the boundaries S m 
of the bodies T> m are C 2 — smooth. 

The heat equation can be stated as follows: 



u t = V 2 u + f(x) in K 3 \ (J V m := w| t=0 = 0, (1.2) 

m— 1 

ii 7v = Q m u on S m , 1 < m < M. (1-3) 



Here N is the outer unit normal to S := II S m , Cm = - m \ x m € 2? m , 1 < m < M, where h(x) is a 

m—l 

continuous function in T>. 
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The aim of this paper is to establish a method for solving both theoretically and numerically the 
many-body heat transfer problem, and to derive an equation for the effective (self-consistent) field in the 
medium in which many small bodies (particles) are embedded. 

There is a large literature on homogenization theory which deals with similar problems, see [lj]-[I| and 
the literature therein. In most cases it is assumed in the homogenization literature that the coefficients of 
the equations are periodic, and there is a small parameter in these coefficients, that the related problems 
in a "periodic cell" have discrete spectrum and are selfadjoint, that the elliptic estimates, such as Rellich 
and Poincare inequalities hold, etc. Our method, which was used in [j|-[l9j. differs in many respects from 
the homogenization theory methods developed in the literature: periodicity assumption is not used, the 
spectrum of the related elliptic problems that we use is continuous, the problems may be non-selfadjoint, 
and our justification of the limiting equation for the effective field is based on a new technique, namely, 
on convergence results for collocation methods (see, Q). The published homogenization techniques are 
not directly applicable to our problem because of the lack of periodicity. 

The main results of this paper are: 

1) Derivation of the linear algebraic system (25) for finding the Laplace transform of the solution to 
problem (2)-(3), 

2) Derivation of the equation (29) for the limiting effective field in the medium when the size of the 
small bodies tends to zero while the total number of these bodies tends to infinity according to the 
distribution law (1), 

3) Derivation of the formula (31) for the average temperature in the limiting medium. 

In Section 2 a derivation of these results is given, and in Section 3 proofs of some lemmas are given. 

2. Derivation of the equation for the limiting effective field 

Denote 

/•OO 

U :=U(x,X) = / e~ xt u(x,t)dt. 
Jo 

Then, (Tj~J) - (fOJ) imply 

- \> 2 U + XU = X^ 1 f(x) in Q, (2.1) 
U N = ( m U on S m , 1 < m < M. (2.2) 



Let 



e --/X\x-y\ ^ 

9{x,y) ■= g(x,y, A) = — -, - / g(x,y)f(y)dy = F(x,X). (2.3) 



4tt\x — y\ X 
Look for the solution to (|2TTj) - (j2T2"j) of the form 

M „ 

U(x, A) = F(x, A) + V / g(x, s)a m (s)ds, U(x, A) := U{x) := U, (2.4) 
m=1 Js m 

where a m are unknown and should be found from the boundary conditions (|2.2[) . Equation (|2.1[) is 
satisfied by U of the form (|2.4[) for any a m . To satisfy (|2.2[) one has to solve equation 

+ A m a m -a m ^ Q m T m a' m = on <S m , 1 < m < M. (2.5) 

Here 



s„ 



U e {x) := Ue t m(x) ■= U{x) - I g(x,s)a m (s)ds, (2.6) 

j(s,s' 
ONs 



dg(s, s') 

T m a m = I g(s,s')a m (s')ds', A m a m = 2 I — — - — a m {s')ds' , (2.7) 



s„ 
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and the known formula for the outer limiting value on S m of the normal derivative of a simple layer 
potential was used. We now apply the ideas and methods for solving many-body scattering problems 
developed in [j| - 0. 

Let us call U e ^ m the effective (self-consistent) value of U, acting on m-th body. As a — > 0, the dependence 
on m disappears, since 

g(x, s)a m (s)ds — > as a — > 0. 



One has 

M 



U(x, A) = F(x, A) + ^ d( x > x m)Q m + J2, x m € T> m , (2.8) 

m— 1 

where 

Qm ■ = / &m{s)ds, 
JS m 

and 

M M 

J 2 := ^2 / [g(x,s') - g(x,x m )]a m (s')ds' , := ^ g(x,x m )Q m . (2.9) 
\Ji\ « \J\\ as a ->• (2.10) 



,,1=1 ^ 5 

We prove (in Section 2) that 
provided that 



Um-^=0. (2.11) 



If (|2.10p holds, then problem (|2. 1[) - (|2.2p is solved asymptotically by the formula 



M 



g(x, x m )Q m , a ->• 0, (2.12) 

m— 1 

provided that asymptotic formulas for Q m , as a — > 0, are found. To find formulas for Q m , let us integrate 
T5"|) over 5 m and estimate the order of the terms in the resulting equation as a — > 0. We get 

/ V 2 U e dx = 0{a 3 ). (2.13) 



Here we assumed that | V 2 W e | = 0(1), a — >• 0. This assumption will be justified in Section 2. 

£ A m° m -a m ds = _ Qm[l + o{l)l a ^ Q (214) 

This relation is justified in Section 2. Furthermore, 

-Uj s U e ds = -Q m \S m \U e {x m ) = 0{a 2 -*), a^0, (2.15) 

where 

\S m \ =0(a 2 ) 

is the surface area of S m . Finally, 

-Cm / ds g(s,s')a m (s')ds' = -C m / ds'o m (.s') / dsg(s,s') 
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= Q m O(a 1_K ), a^O. (2.16) 
Thus, the main term of the asymptotics of Q m is 

Qm = ~Cm\S m \lA e {x m ). (2-17) 

Formulas ((2~lTj) and (|2~T2"j) yield 

M 

U{x,X)=F(x,X)-J2 ( m \S m \U e {x m ,X), (2.18) 

m— 1 

and 

M 

U e (x mi A) = F(x m , A) - ^2 )Cm'|<Sm'|^e(a;m' 1 A). (2.19) 

m'^m,m' — 1 

Denote 

I4 e (x m , A) := i/ m , 

F (%m i A) • F rn: ^(^mi ^m') ■ — 9 mm' i 

and write (|2 . 19[) as a linear algebraic system 

W m = F m - a 2 ~ K 9m m 'h m >c m ,U m ,, 1 < m < M, (2.20) 

where ^ 

Consider a partition of the bounded domain T>, in which the small bodies are distributed, into a union 
of P << M small nonintersecting cubes A p , 1 < p < P, of side b >> d, b ~ b(a) — > as a — >• 0. Let 
x p G Z\ p , |Z\ p | = volume of A p . One has 

m p 

CL ^ ^ Qmm'hm' C-m'Hm' ^ ^ ^ Qpp' hp' CpflApf ^ ^ 1 

m' — l^m'^m p' — l^p'^p x m /£A p i 

= Y,9pp>hp>c p ,U p ,N(x p ,)\A p ,\[l + o{l)}, a^O. (2.21) 
p'^p 

Thus, (j2~20)) yields 

p 

U v = F p - g pp ,h p ,c p ,N p ,U p ,\A p ,\, 1< P <P (2.22) 

We have assumed that 

hm> = V[l + »(l)]i c m / = v[l + o(l)], W m / =W P 4l + o(l)], a^O, (2.23) 
for x m < G Zip'. This assumption is justified if the functions h(x) 1 U(x, A), 

/ \ i • I ^m' I 

c(x) = hm — — , 

and 7V(:r) are continuous. The function h(x) and N(x) are continuous by the assumption. The continuity 
of the U(x, A) is proved in Section 2, and the continuity of c(x) is assumed. If all the small bodies are 

35 



621" — 2012/7/4 — 0:14 — page 36 — #5 



A. G. Ramm Heat transfer in a medium 

identical, then c(x) = c = const. 

The sum in the right-hand side of (|2.22l) is the Ricmannian sum for the integral 

p r 

N{x p ,)U p ,= / g(x,y)h(y)c(y)N(y)U(y,X)dy (2.24) 

p>=l >p >jip 

Therefore, linear algebraic system (|2.22p is a collocation method for solving integral equation 

U{x, X) = F(x, A) - / g{x, y)h{y)c{y)N{y)U{ Vl X)dy. (2.25) 
Jv 

Convergence of this method for solving equations with weakly singular kernels is proved in Q . 
Applying the operator —V 2 + A and then taking the inverse Laplace transform of (|2.25p yields 

u t = Au + f{x) - q{x)u, q(x) := h(x)c{x)N(x). (2.26) 

One concludes that the limiting equation for the temperature contains the term q{x)u. Thus, the em- 
bedding of many small particles creates a distribution of source and sink terms in the medium, the 
distribution of which is described by the term q(x)u. 

If one solves equation (|2.25p for U(x, A), or linear algebraic system (|2.22[) for U P (X), then one can Laplace- 
invert U(x, A) for U{x, t). Numerical methods for Laplace inversion from the real axis are discussed in 0] 

-M- 

If one is interested only in the average temperature, one can use the relation 

i r T 

lim — / u(x, t)dt = lim XU{x, X) := ip(x). (2.27) 

T->oo 1 Jq \~>0 

Relation f|2 . 2T[) is proved in Section 2, which holds if the limit on one of its sides exists. The limit on the 
right-hand side of (|2.27p can be calculated by the formula 

i>{x) = (I+ B)~ V, if = [ -J(y)dy. (2.28) 

J 47r|x - y\ 



Here, B is the operator 



I #W 4 q{x) ._ h{x)c{x)N{x y 



From the physical point of view the function h(x) is non-positive because the flux — Vu of the heat flow 
is proportional to the temperature u and is directed along the outer normal ./V: —un = h\u, where 
hi = —h > 0. Thus, q < 0. It is proved in [ll[ - [l2[ that zero is not an eigenvalue of the operator 

— V 2 + q(x) provided that q(x) > and q = 0( ) as \x\ — > oo, e > 0. In our case, q(x) = outside 

|x| /+e 

T> 1 so the operator (I + B)~ l exists and is bounded in C(T>). Let us formulate the basic result we have 
proved. 

Theorem 1. Assume hl.l}) , \2.11\) , and h < 0. Then, there exists the limit U(x,X) ofU e (x,X) as 
a — > 0, lA{x 1 X) solves equation V2.25)) , and there exists the limit \2.2ty , where ip(x) is given by formula 

WW- 

3. Proofs of some lemmas 

Lemma 3.1. Assume S2.ll)) . Then relation 12. 1 0|) holds. 
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Proof. One has 



J hm := \g(x,x m )Q\ = 0(' l%lf — p. ) < \ X - Xm \ > d . 

4n\x-x m \ J \x - x m \ 



-v\\x— x m \ 



Is 47r|a; — x 
where \x — x rn \ > 2. and the inequality 



max 



VXa, 



\x x rn \ 



\a m {s')\ds' < O 



\Qm\a 



A>0 \X — X r 



was used. The \Q m \ ^ 0. In fact, <r m keeps sign on S m , as follows from equation (|2.5[) as a — > 0. 
It follows from - (JO) that 



^2, 



J: 



l,m 



< o 



From (|3.3[) by the arguments similar to the given in [13( one obtains ([2~10| . 

Lemma 3.2. Relation {2.1$ holds. 

Proof. Let us justify relation (|2.14j) . As a — >• 0, one has 

9 e-^l 8 " 8 '! 9 1 



d e -^\ s ~ s '\ — i 
dN~ s 4tt|s - s'| = 9A^4tt|s- s'| + 9iV s 4tt|s-s'| 



It is known (see Q) that 



ds 



d 1 



/s m is m dN s 47r|s - s' 
On the other hand, as a — > 0, one has 

L ds L am{s ' )ds ' 

The relations ([33]) and justify ([2TT4]) . 
Lemma 3.3. Relation $2.27]) holds. 
Proof. Denote 



-a m (s')ds' = - I a m {s')ds' = -Q, 



< \Qr 



ds- 



1 _ e -V\\s-s'\ 



4tt\s - 



o{Q m )- 



1 f* 

- / u(t)dt := v(t), u(a) := 



Then 



e(A) 



u(a) 
a CT 



f/rr 



by the properties of the Laplace transform. Assume that the limit v(oo) := Woo exists: 

lim v(t) = Voo. 

t— >oo 

Then, 

/>oo 

Uoo = lim A / e~ xt v(t)dt = lim Xv(X). 

A-s-0 '- 



A->0 
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(3.3) 
□ 



(3.4) 
(3.5) 

(3.6) 
□ 



(3.7) 
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Indeed, 



so 



A / e~ xt dt = 1 



lim A / e' M (v{t) - Voo )dt = 0, 



and (|3.7|) is verified. One has 



Let us check this: 



[°° A 

lim A"D(A) = lim / —u(a)da = lim Au(A). (3-8) 
A-s-0 A-s-0 J \ a A-s-0 

/*°° A f°° A 

lim / —u((j)da = lim / -^-<ju(a)d<T = lim <ju(<t), (3.9) 
A-s-o J x a A-s-o J x a 2 cr^o 



where we have used the relation / — da = 1. 



A 



a 2 



Alternatively, let a =7. Then 



A 1 f 1/x 1_1 1 f u 1 _ 1 

-^r(ju((j)d(j — — — I —u( — )d-f=— / —u(—)d'y. (3.10) 

a CT V A Jo 7 7 W 7 7 

If A — > 0, then cj = A -1 — > 00, and if i/» := 7~ 1 m(7~ 1 ), then 

lim — / ipdr/ = ip(oo) = lim 7 _1 u(7 _1 ) = lim au(a). (3-11) 

uj— >oo Ul / n 7— >oo CT— ^0 



□ 
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